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ABSTRACT 

We study the dependence of the A<BH-A/host relation on the redshift up to z = 3 for 
a sample of 96 quasars the host galaxy luminosities of which are known. Black hole 
masses were estimated assuming virial equilibrium in the broad line regions (Paper I) , 
while the host galaxy masses were inferred from their luminosities. With this data we 
are able to pin down the redshift dependence of the AiBH-Afhost relation along 85 per 
cent of the Universe age. We show that, in the sampled redshift range, the A^BH-^host 
relation remains nearly unchanged. Once we take into account the aging of the stellar 
population, we find that the AfBH/A^host ratio (T) increases by a factor ~ 7 from 
z = to z = 3. We show that T evolves with z regardless of the radio loudness and of 
the quasar luminosity. We propose that most massive black holes, living their quasar 
phase at high- redshift, become extremely rare objects in host galaxies of similar mass 
in the Local Universe. 
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1 INTRODUCTION 

Many pieces of evidence suggest that supermassive black 
holes (BHs) and their host galaxies share a joint evolution 
throughout Cosmic Time. In part icular: i- the evolution of 
the quasar luminosity function (iDunlop fc Peacockl Il990l ; 



iFontanot et all 120071 : ICroom et al 
the trend of the star formation 



2009) 



closely matches 
density t hrough Cos- 
mic ages ( Ma dau Pozzetti fc Dickinson! [l998h : ii- massive 
black holes are found in virtually all massive galaxies 
(|Kormendv fc Richstonelll995l ; lDecarli et alfeOOTT l; in- their 
mass (AIbh) is tightly correlated with the large scale 
properties (stellar velocity dispersion, a* ; l uminosity, Z/host ; 
mass, At host) of th eir host galaxies (see iFerraresg 120061 : 
iGultekin et al. 2009, for recent reviews on this topic). 

When and how these relations set in, and which are 
the physical processes responsible of their onset are still 
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from a theoretical (e.g., Silk & Rces 


199ci_ Kind 20051: 


Wvithe & Loebl 20061: Robertson et al. 


2006; iHoDkins et alj 


20071: iMalbon et al. 2007) and an 


observational point of view 


(McLure et al. 


2006:|Pene: et all l2006alfbl; ISalviander et all 
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Probing BH-host galaxy relations at high redshift is 
extremely challenging. Direct measurements of Mbh from 
gaseous or stellar dynamics require observations capable to 
resolve the BH sphere of influence, which are feasible only 
for a limited number of nearby galaxies. The only way to 
measure Mbh in distant (> few tens of Mpc) galaxies is 
to focus on Type-1 AGN, where Mbh can be inferred from 
the width of emission lines Doppler-broadened by the BH 
potential well and from the AGN continuum luminosity 
(see, e.g., IVestergaardl 120021 ). assuming the virial equilib- 
rium. Quasars represent therefore the best tool to probe 
A4bh at high redshift, thanks to their huge luminosity. In- 
deed, large-field spectroscopic surveys allowed the estimate 
of .Mbh in several thousands of obje cts (see, for instance, 
IShen et al.ll2008al ; lLabita et al.ll2009 a b). On the other hand, 
the AGN light in quasars outshines the emission from the 
host galaxies, making their detection challenging. Only re- 
cently limitations due to intrinsic (e.g. the Nucleus-to-Host 
galaxy luminosity ratio, N/H) and extrinsic (e.g., the angu- 
lar size of the host with respect to the angular resolution of 
the observations) effects could be overcome. Optical images 
from the Hubble Space Telescope (HST) and NIR observa- 
tions both from space and ground-based telescopes could 
resolve ~ 300 quasar host galaxies up to z ~ 3 (see, e.g., 
iKotilainen et al.1l2009l . and references therein). Preliminary 
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studies suggest that, for a given galaxy mass, black holes in 
high- z AGN are more ma s sive than their low -z counter-parts 
(e.g., iMcLure et alll2006l : IPeng et al.ll2006al lbh. 

A number of limitations potentially affect the studies 
of the A^BH-Alhoat relation through Cosmic Time: 1) All 
the works to date use different .Mbh proxies as a function 
of redshift (i.e., based on different broad emission lines: usu- 
ally UP at z < 0.5, Mgn for 0.5 < z < 2 and Civ at 
z 1.6). 2) Selection biases related to luminosity or flux 
limits, to the sampled N/H, to the steepness of the bright 
end of the galaxy and quasar luminosity functions may hin- 
der the study of the evolution of t he -MBH-.Mhost relation 
(see, for instance. lLauer et al.ll2007r ). 3) As the properties of 
quasar host galaxies are directly observed only in a limited 
number of objects, poor statistics usually affect the available 
datasets. 

This study represents a significant effort in overcoming 
all these limitations: Thanks to UV and optical spectra o f 
low-redshift quasars l|Labita et al.ll2006l : iDecarli et alj|2008h 
and to optical spect ra of mid- and high-redshift quasars 
l|Decarli et al.ll2009bl . hereafter Paper I) we can probe .Mbh 
using both high and low ionization lines in a wide range of 
redshifts, for the largest dataset adopted so far in this kind 
of studies. 

In Paper I we presented the sample and we inferred BH 
masses. Here we describe the data sources for the host galaxy 
luminosities, we infer .Mhost and we address the evolution 
of the .Mbh -.Mhost relation. 

Throughout the paper, we adopt a concordance cosmol- 
ogy with H = 70 km/s/Mpc, Q m = 0.3, fi A = 0.7. We con- 
verted the results of other authors to this cosmology when 
adopting their relations and data. 



2 L host AND .M host IN DISTANT QUASARS 

For the purposes of this work, we have to define a homoge- 
neous compilation of quasar host galaxy luminosities from 
data available in the literature. After that, we use the rest- 
frame _R-band lumin osities to infer the hos t galaxy stellar 
masses. We refer to iKotilainen et all (2009) for a detailed 
discussion of technicalities in the luminosity estimate of 
quasar host galaxies from high-resolution imaging. 



2.1 Host galaxy luminosities from the literature 

For a complete list of data sources, we refer to the sample 
description in Paper I. Apparent magnitudes in the filters 
of the observations are converted to rest-frame 7?-band ab- 
solute magnitude as follows: 



Mb 



m f -5logD h (z)-C{z)-A f 



(1) 



where / is the original filter of the observations, -Dl(z) is 
the luminosity distance of the quasar in the cosmological 
frame we adopted, C(z) is a term accounting for filter and k- 
correc tion, as derived by ass uming an elliptical galaxy tem- 
plate (|Mannucci et al.ll200ll ). and Af is a term accounting 
fo r the Galactic extinction, as derive d from the Hi maps 
m ISchleeel Finkbeiner fc Davis! (|l998T ). We remark that, in 
order to minimize filter and colour corrections, we selected 
observations performed using filters roughly sampling the 
rest-frame i?-band. Moreover, the i?-band luminosity is only 



marginally sensitive to the age of the stellar content. Thus, 
uncertainties in C(z) due to the chosen host galaxy template 
are negligible (< 0.1 mag) for the purposes of this work. 

Low-z data taken with the HST-Wide Field Camera 
have been analyzed by many authors, and different m/ esti- 
mates are available for the same ob ject and on t h e sam e 
data. In partic u lar, t h e studies of Bahcall et al.l l|l997h ; 
lHamilton et all (|2002r ): iDunlop et all (|2003h si gnificantly 
overlap onto the recent re-analysis presented by Kim et all 
(|2008al lbh. When comparing the reported apparent host 
galaxy and nuclear magnitudes, we find that the average off- 
set is usually negligible (^$ 0.2 mag), but a significant scatter 
is present (rms~ 0.3 — 0.5 mag). When more than one esti- 
mate of rrif was available, we adopted the most recent one. 
No images of the mid- and high-z quasars in our sample 
were analyzed independently by different groups, thus no 
superposition happens for these objects. 



2.2 Host galaxy masses 

In order to infer the stellar mass from the host galaxy lu- 
minosity, we have to adopt a stellar _R-band Mass-to-Light 
ratio and consider its dependence on Cosmic Time. If the 
majority of the stellar population of massive galaxies did 
form at high redshift, as suggested by several pieces of 
evidence dGavazzi et al.ll2002l: [Thomas et al.ll2005l iRenzinil 
l2006l : ICirasuolo et al.ll2008l : ICappellari et al.ll2009h , one may 
assume that the Mass-to-Light ratio passively evolves from 
the formation (z = Zburst) to the present age. On the 
other hand, if quasar host galaxies suffer intense star for- 
mation episodes from z = 3 to z = (for instance, due to 
merger events), the evolution of the stellar Mass-to-Light 
ratio becomes more complex and is, in principle, differ- 
ent fro m object to object. Fo r the sake of simplicity, fol- 
lowing IKotilainen et al] (|2009h we will consider here only 
the scenario of a passively evolving stellar population with 
Zburst = 5. This is justified by the selection of quasars with 
bulge-dominated host galaxies, where old stellar populations 
are expected. Furthermore, as we will discuss in section 3J 
this assumption is conservative with respect to the main re- 
sults of our study. 

With this caveat in mind, we find that the redshift de- 
pendence of the host galaxy luminosity observed, e.g., in 
IKotilainen et al] (|2009| ) is practically removed when we take 
into account the evolution of the stellar population. The stel- 
lar mass of the host galaxies in our sample is nearly constant, 
with an average value of few times 10 1 M©. 

Table Q] lists our final estimates of the host galaxy lu- 
minosities and masses for the quasars in our sample. 



3 EVOLUTION OF THE M B u-L host , M hos t 
RELATIONS 

In Figure [T] we compare our .Mbh e stimates with the pre- 
dictions of the iBettoni et all (2003) relation, defined on 
z ~ galaxies, and with the expectations in the case of 
a fixed -MBH/-Mhost= 0.002 ratio, as observed in Local inac- 
tive galaxies (see, e.g.. iMarconi fc Huntll2003l ). 

The A^BH-ihost relation appears rather insensitive to 
the Cosmic Time, independently on which line is adopted 
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Table 1. The sample in this study. (1) Quasar name. (2) Catalogue redshift. (3) Host galaxy ij-band absolute magnitude (not corrected 
for st ellar aging). (4 ) Host gala xy stellar mas s . (5) R efer ence for the host galax y lu minosity: 1- low-z H ST-b ased observation ( see P aper 
II: 2- IFalomo et all d2004l): 3- iKukula et all (l200lh : 4- iRidgwav et alj <200ll) : 5 - IFalomo et alJ ll2005ll: 6- IKotilainen et al.l d2009l): 7- 
iKotilainen et alj J2007I' ): 8- lHvvonen et al.l d2007al) ; 9- [Kotilainen et al.l lll998l) ; 10- IKotilainen et alj ||2000|) ; 11- lDecarli Treves fc Falomol 
ll2009ah . 
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in the virial estimate of Mbh- When correcting for the evo- .Mhost and their ratio F are plotted all together as a function 
lution of the stellar Mass-to-Light ratio, we find a clear in- of redshift. The linear best fit of log V is: 
crease (~ 0.7 dex) of the A^BH/A^host ratio with respect 
to what observed in the Local Universe. In Figure [2] Mbh, 



logT = (0.28 ± 0.06) z - (2.91 ± 0.06) 



(2) 
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log M host [M ] 

Figure 1. The Mbh~ ^host an d A / lBH~A^host relations in three different redsh ift bins. Squares (tri angles, circles) mark quasars in 
which Mbh is derived from H/3 (Mgn, Civ). The reference (solid) line is the iBettoni et al. (2003) relation (upper panels) or the 
-A / 'BH/A / fhost = 0.002 case (lower panels). The dotted line is the best fit to the data, assuming the same slope of the rest-frame relations. 
No significant redshift evolution is observed when comparing Mbh with the observed host galaxy luminosities. On the other hand, a 
clear offset is apparent in the .Mbh --Mhost relationship as a function of the redshift. 



Table 2. Slope of the best linear fit of log r = az + /3, for the 
whole sample and various subsamples (column 1). The number 
of objects in each subsample in given in col. 2. Uncertainties on 
the slope (col. 3) are analytically derived from the Least Square 
minimization criterion. 



Subset 


N.obj. 


a 


(1) 


(2) 


(3) 


All 


96 


(0.28 ±0.06) 


RLQs 


48 


(0.24 ±0.11) 


RQQs 


48 


(0.31 ± 0.07) 


M * > M Y > M* - 1 


33 


(0.32 ± 0.10) 


-26 > My > -27 


31 


(0.32 ± 0.11) 


N/H>5 


43 


(0.33 ± 0.11) 


N/H<5 


53 


(0.25 ± 0.07) 



suggesting that galaxies with similar stellar masses harbour 
BHs ~ 7 times more massive at z = 3 than galaxies at z = 0. 

In Figure [3] we study separately Radio Loud (RLQs) 
and Radio Quiet quasars (RQQs) , finding that both samples 
are consistent with the log T—z relation found for the whole 



sample. The only remarkable difference is in the offset, in 
the sense that, at any redshift, both black holes and host 
galax ies in RLQs are ~ . 2 dex more massive than in RQQs 
(e.g. , iDunlop et ail 120031 ; lLabita et all l2009d ) . Table [2] and 
Figure U report the slopes of the best linear fit of log T as a 
function of z in each subsample. 

3.1 Is the trend of V an artifact? 

In this section, we discuss the possible effects that could 
bias the estimate of F, in order to probe the reliability of 
the trend observed in Figures [TH3 

3.1.1 The luminosity function bias 

lLauer et al.l (|2007l ) showed that, because of the steepness of 
the bright end of the galaxy luminosity (mass) function and 
the presence of intrinsic scatter in the A^BH-ihost (A^host) 
relation, very massive BHs are preferentially found in rela- 
tively faint (less massive) galaxies rather than in extremely 
bright (massive) galaxies, which are very rare. Since high-z 
samples are dominated by massive objects, the bias increases 
with z, possibly mimicking an evolution in F. 
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-4 - 

12 3 

Redshift 

Figure 2. The redshift dependence of AIbh (top panel), Afhost 
(middle panel) and their ratio T (bottom panel). The symbol code 
follows Figure [l] The best linear fits are plotted. The average 
points with rms as error bars of the H/3 subsample (big square), 
of the low- and high-z C IV data (big circles) and of the Mg II data 
with redshift < 1 and > 1 (big triangles) are also shown. 



In order to quantify the relevance of this bias, we assume 
that Mbh mostly depends on the quasar luminosity. This 
is consistent with the relatively small range of Eddington 
ratios we sample. If cr M , the cosmic scatter of the A^BH-khost , 
is constant in Lhost> at a given redshift the bias depends on 
the shape of th e luminosity funct ion of quasars, ^t(M) (see 
equation 25 in lLauer et aUl2007h . We assume the quasar 
luminosity f unction and its pu rely-luminosity evolution as 
reported bv lBovle et all (|200Ch , basing on the 2QZ survey: 
M*(z) = -22. - 2. 5 (1.34 2-0. 27 z 2 ). As a consequence, as 
long as we sample the same range of the quasar luminosity 
function at any redshift, the bias is kept constant. A constant 
bias is irrelevant for the purposes of our study, since our 
main aim is to probe the redshift dependence of the BH~ 
host galaxy relations, not their absolute normalization. 

From Figure 1 of Paper I, we note that the constant 
luminosity cut at -26 > My > -27 and the M* > M v > 
M* — 1 cut roughly braket the objects in our sample over 5 
magnitudes in My. Hence, if we consider the whole sample, 
the bias on A^bh will lie within the expectations from these 
two cases. Our estimate of the redshift evolution of the bias 
in these two cases is plotted in Figure E for two diff erent 
values of <t m , namely 0.39 from iBettoni et al.l (|2003l ) and 



-l F 




12 3 

Redshift 



Figure 3. The redshift dependence of Y for Radio Loud (top 
panel) and Radio Quiet quasars (bottom panel) separately. The 
symbol code is the same as in Figure[2] The number of objects in 
each subsample is also provided in parenthesis. 



the conservative value <r M = 0.5. We conclude that the bias 
accounts for < 0.11 dex « a factor 1.3 moving from z — to 
z = 3. As the observed dependence of F is ~ 6 times larger, 
it cannot be explained in terms of this selection effect. As a 
further check, Figure [6] shows the log T—z plane only for the 
objects lying in the two luminosity cuts considered in this 
discussion. The observed trend is unchanged, independently 
on the adopted luminosity cut (see Table [2] and Figure [4}. 

3.1.2 The effects of the N/H ratio 

All the objects in our reference sample are selected on the 
basis of their total luminosity, which is dominated by the 
nuclear light in quasars. This possibly introduces a bias in 
the sense that the higher is the Nuclear-to-Host luminosity 
ratio (N/H), the harder is the measure of the host galaxy 
luminosity, especially at high-z. In Figure [7] we show that 
the redshift dependence of F in objects with high- and low- 
N/H is similar. Moreover, we stress that if we include the 
unresolved quasars to this analysis, the trend would be even 
steeper, as they all lie at high-z. Analogously, this argu- 
ment can be applied for possible contaminations from disc- 
dominated galaxies at high redshift, where the morphology 
classification may be more doubtful. In this AIbh is 

sensitive to the bulge mass rather than to the total galaxy 
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— -26>M v >-27, ct m = 0.5 

— M,>M V >M,-1, ct m = 0.5 
26>M v >-27, a = 0.39 
M,>M„>M,-1, a = 0.39 




Figure 4. The slope of the log T versus 2 linear fit in our data, 
for each subsample. Error bars are the 1-cr uncertainties as de- 
rived from the fit algorithm (see Table [2} . All the subsets have 
consistent slopes around ~ 0.3 dex. A non-evolving scenario (hor- 
izontal, dotted line) mismatches with the observations in all the 
cases. 



mass, we should consider smaller values of A^host for such 
galaxies, which would increase the value of F at high z. 



Figure 5. The bias on the prediction of .Mbh from the -Mbh~ 
Thost relation with respect to the expectation from the luminos- 
ity functions of galaxies and quasars, plotted as a function of 
redshift. The bias estimates are obtained by integrating the lu- 
minosity function of quasars over the adopted luminosity cuts: 
-26 > My > -27 (dot-dashed and dotted lines) and M, > 
My > M* — 1 (dashed and solid lines). We plot the limit cases 
with <t m = 0.5 (dot-dashed and dashed lines) and = 0.3 (dot- 
ted and solid lines). We note that, in the worst case, the bias 
increases of 0.22 dex (that is, a factor 1.66) from z = 2.5 to 

2 = 0. 



3.1.3 The role of radiation pressure 

iMarconi et al 1 (|2008l . I 20091 ) suggest that the virial estimates 
of .Mbh may yield lower limits to the true BH mass, as the 
radiation pressure is not taken into account. As long as the 
BLR clouds are virialized, a correction can be applied by 
adding a term depending on the BLR column density ATr 
and the quasar luminosity. There is still no strong constrain 
on the values of Nh- X-ray variability studies performed 
on nearby, lower luminosity AGN su ggest that the column 
densities may be re latively high (e.g. iTurner fc Millerll2009l ; 
iRisaliti et al.ll2009l . and references therein), thus preventing 
radiation pressure from sustaining BLR clouds motion. How- 
ever, as a clear comprehension of the radiation pressure role 
in the BLR is still missing, especially in the most luminous 
AGN, we limit our discussion to the following consideration: 
Since the average luminosity of our data increase with z, the 
radiation pressure effect is expected to become more severe 
at higher z, leading to an even steeper trend of F than the 
one observed in Figures [T] and [5] 



4 DISCUSSION 

4.1 Comparison with previous results 

As a key result of this analysis, we find that the -MBH/.Mhost 
ratio significantly increases with redshift. Hereafter we com- 



pare these findings with those of other studies available in 
the literature. 

iMcLure et ail (|2006n match the average trend of .Mbh 
observed in 38 Radio Loud quasars at z < 2 with the typical 
stellar masses of massive radio galaxies in the same redshift 
bins. This approach relies on the assumption that quasar 
host galaxies are comparable, at any redshift, with massive 
radio galaxies. The major caveat here is that their results 
may be biased by the different histories of quasar and radio 
galaxies (for instance, note that the luminosity functions of 
AGN evolves differently for vario us luminosity subclasses). 
Nevertheless. lMcLure et"al1 (|2006l ) find an increase of V com- 
parable to the one observed in the present study (see Figure 
[SJ. Our results extend these findings to RQQs and beyond 
the peak age of q uasar a ctivity. 

IPeng et alj (|2006al lbl) address the evolution of the 
-MBH--Mhost relation as a function of redshift in ~ 20 low-z 
quasars and in ~ 30 high-z lensed quasars imaged with the 
HST. Their data show a larger scatter than ours, possibly 
due to uncertainties in the modelling of the lens mass distri- 
bution and the lens light subtraction. They find that there is 
practically no evolution in the A^BH-^host relation. On the 
other hand, when correcting for the evolution of the stellar 
population, they find an excess (a factor 3-6) in the .Mbh 
values at high-z with respect to the prediction from the lo- 
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Figure 6. The same as Figure [3] for objects with M* > My > 
M* — 1 (top panel) and with —26 > My > —27 (bottom panel). 
No significant difference is observed in the two subsets. 



Figure 7. The same as Figure [3] for objects with N/H>5 (top 
panel) and N/H<5 (bottom panel). Again, no significant differ- 
ence is reported. 



cal A^BH-A^host relation, in qualitative agreement with our 
find ings. 

iMerloni et al.l ([2009) study the AtBH-A'lhost ratio in a 
sample of 89 type-1 AGN with 1 < z < 2.2 from the zCOS- 
MOS survey. Black hole masses are derived through the stan- 
dard virial assumption, while the host galaxy luminosities 
and stellar masses are inferred from multi-wavelength fitting 
of the spectral energy distribution of the targets (with no di- 
rect information about the morphology of the galaxies). This 
technique is effective with intermediate to low-luminosity 
AGN, while it cannot be applied to quasars as bright as ours, 
where the nuclear light overwhelms the galaxy contribution. 
They find that the average F is higher than what observed in 
the Local Universe, the excess scaling as (1 + z) ' 74 , consis- 
tently wi th the trend observe d in our data in the same red- 
shift bin. IJahnke et all (2009) observed 10 of the targets in 
Merloni's sample with the HST and independently derived 
host galaxies luminosities with a procedu re similar to the 
one a dopted in our data sources (see, e.g.. iKotilainen et al.l 
2009). They find no evolution in the Mbh~M host (total) ra- 
tio. However, clues of the occurrence of discs are present, 
thus the MBH-Mhost (bulge) ratio is expected to evolve as 
(1 + z) 12 , in agreement w ith our findings. 

iBennert et all l|2009h address the A'fBH-^host relation 
in a sample of 23 Seyfert galaxies with 0.3 < z < 0.6. They 
study the morphology of the host galaxies of these objects 



using NIR observations from the HST. A careful modelling 
is adopted in order to disentangle nuclei, bulges and disc 
components. Black hole masses are derived in a way similar 
to that presented in paper I. Once corrected for the evolution 
of the stellar population, they find T oc (1 + z)' 14±0,2 ' ) , in 
good agreement with our results. 

Additional indication of an evolution of F comes from 
the relation of A(bh with the stellar velocity dispersion, 
a,, of the host galaxy. In particular, it is remarkable that 
these works suggest that the higher is the redshift, the 
more massive is t he bla ck hole for a given cr*. For instance, 
ISalviander et all l|2007h use the width of the [O in] narrow 
emission line as a proxy of the stellar velocity dispersion and 
study the AIbh-o"* relation in a sample of ~ 1600 quasars 
up to z — 1.2 taken from the SDSS. They find that A(bh at 
high redshift are ~ 0.2 dex larger than what expected from 
local A(bh-o"* relation. A smaller evolution (< 0.1 dex), al- 
beit wi th small significance, is also proposed by IShen et al.l 
i|2008br i, based on a sample of 900 Type-1 AGN with z < 0.4. 
More recently, IWoo et all l|2008h and Woo et al. (in prepa- 
ration) address the A^bh-ct* relation in Seyfert galaxies up 
to z ~ 0.6, and find an overall AIbh excess at high red- 
shift with respect to the prediction from low-z relationships. 
These findings support our results, notwithstanding the dif- 
ferent characteristic luminosities, morphologies and stellar 
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Figure 8. The average values of T for the quasars in our sample 
(filled symbols), together with the linear best fit (dotted line). 
Error bars are the 2-cr uncertainties in the average values for each 
data bin. For a comparison, the tren d found in the study of radio 
loud AGN by iMcLure et"aD (|2006h is plotted as a dashed line. 
We also plot the average T values of the 51 lensed and non-lensed 



quasars fr omlPe ng et al 



AGN from I Merloni et al 



( 2006a. h, empty circles) and of the 89 
(2009, empty triangles). The horizontal 
solid line represents the constant T = 0.002 case. All together the 
data depict a clear increase of T with z. 



contents of the sampled targets with respect to those exam- 
ined in our analysis. 

It is interesting to note that the z = 6.42 quasar 
SPSS J1148+5251 has a black h ole mass of few x 10 9 M 
(|Willott. McLure fc Jarvisl 120031) and a dynamical mass o f 
the host galaxy of ~ 5 x 10 10 M Q (| Walter et al.ll2003l , l2004T) , 
yielding F ~ 0.1, which is in agreement with the extrap- 
olation of our results at that redshift (r « 0.13) and well 
beyond the V — 0.002 value observed in the Local Universe. 

These results as a whole support a picture where, for 
a given quasar host galaxy, its central black hole at high 
redshift is 'over-massive' with respect to its low- z counter- 
parts. This picture is also consistent with the constrains on 
the A^BH-A^host evolution derived from the comparison be- 
tween the galax y stellar mass fun ction and the quasar lumi- 
nosity function (|Somervillell2009l ) . 



4.2 Why does T evolve? 

The interpretation of the observed evolution in the A^bh- 
A4host ratio is challenging. Exotic scenarios involving black 
hole ejection from their host galaxies due to gravitational 
wave recoil or to 3-body s catter may be applicable for few 
peculiar targets (e.g., see iKomossa. Zhou fc Lul 120081 . but 
see also Bogdanovic et al., 2009, Heckman et al., 2009 and 
Dotti et al., 2009 for alternative explanations), but they are 
not applicable to the general case. Thus, if high-z quasars 



are destined to move towards the local A^BH-Alhost, the 
unavoidable consequence of our results is that, at a given 
A4bh, galaxy masses increase from z — 3 to the present age. 
Hereafter we sketch three possible basic pictures for that. 
We also present an alternative scenario, in which the fate 
of high-z quasars may be different, the remnants of high-z 
quasars keeping high F values down to the present age. 

Galaxy growth by mergers - A first scenario involves 
substantial mass growth of quasar host galaxies through 
merger events. It is remarkable that strong gravitational 
interactions may trigger intense gas infall in the centre of 
galaxies and may even lead to the activation of BH accre- 
tion. This i s observed in a number of relatively low-redshift 
AGN (e.g.. iBennert et al.ll200ci I2009T ) showing dense close 
environments or disturbed morphologies, and confirmed 
by the presence of young stellar populations in some 
quasar host galaxi es ijjahnke. Kuhlbrodt fc Wisotzkil 120041 : 
IJahnke et al.ll2007l ). Two arguments disfavour this scenario. 
First, theoretical models bas ed on the structure e volution 
in a ACDM cosmology (e.g.. IVolonteri et aJjl20o3 ) predict 
that a massive galaxy experience only few (1-2) major 
merger events from z = 3 to z = 0. However, our study 
shows that a factor ~ 7 increase of the stellar mass of the 
host galaxies is required from z = 3 to z = 0, which means 
that the host galaxies have to suffer ^ 3 major mergers 
in this redshift range. Secondly, several pieces of evidence 
suggest that massive inactive galaxies, as well as quasar 
host galaxies, have already formed/assembled the majority 
of th eir mass in very remo te Cosmic epochs (z > 3; see, 
e.g., iKotilainen et al] 120091 . and references therein). The 
stellar population may experience episodic rejuvenation, 
but this only marginally affects the mean a ge of the stellar 
conten t: Th e stellar shells observe d, e.g., by Canalizo et al] 
(|2007h and IBennert et al] j2()08) in low-redshift quasar 
host galaxies account for 5-10 per cent of the total 
stellar population. Similarl y, in a comparison with inactive 
galaxi es of similar mass, IJahnke. Kuhlbrodt fc Wisotzkil 
(2004h find that quasar host galaxies are on average only 
0.3 mag bluer. If the galaxies enter the quasar phase ~ 1 
Gyr after the activation of the starburst, as suggested by 
the authors of that study, then the involved mass is ~ 30 
per cent of the initial mass of the galaxy. Moreover, if the 
quasar host galaxies contain a significant fraction of young 
stellar populations, then the Mass-to-Light ratio would be 
smaller. Therefore, young host galaxies at high-z would 
yield a F-z relation even steeper than that reported in 
Figured 

Stellar mass growth through gas consumption - Another 
possible interpretation is that high redshift quasar host 
galaxies are gas rich, and form a significant fraction of their 
stellar content in relatively recent Cosmic epochs. This is 
consistent with a picture in which the black hole mass is 
somehow sensitive to the energetic budget of the galaxy 
or its dy namical mass rather than its stellar mass (see for 
instance iHopkins et "all 120071 ) . This scenario is disfavoured 
as all the quasar host galaxies in our sample are massive 
elliptical, and the stellar content of these galaxies is usually 
old. Moreover, if significant star formation occurred in 
quasar host galaxies in the redshift range explored in this 
work, the evolution of V would be much steeper, making 
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this scenario even less realistic. 

Evolution of the fundamental plane - A number of studies 
suggest that inactive, massive elliptical galaxies were more 
compact in the high redshift tha n in the Local Univers e 
(Trujillo et al. 2006; but see also ICappellari et all 120091) . 
In particular, an evolution of the Faber-Jackson relation is 
predicted, implying that the higher is z, the higher is the 
galaxy velocity dispersion a„. If Mbh consta ntly regulates 
the host galaxy <r* (e.g., ISilk fc R ccs 1998), so that the 
-Mbh-o", relation does not evolve significantly, then even 
a small (a factor ~ 1.6) increase of cr, for a given galaxy 
would yield to an excess of a factor 7 in F. The main 
limit of this picture is that studies of the evolution of the 
Mbh~o~, relation through redshift do find an increase of 
the average Mbh for a given a*, when moving from the 
Loca l to high-z Universe l|Salviander et al.ll2007l ; IWoo et al.l 
l200Sl : iBennert et ai1l2009h . 

Remnants of high-z quasars as rare outliers - We propose 
a scenario in which the local counter-parts of high-z quasars 
are high-mass outliers above the A^BH-Alhost relation. 
The more massi ve is the BH, the earlier it exper iences 
its quasar phase (|Marconi et alj|2004l ; lMerloni| [2004'l. Our 
study shows that these objects have higher expected V, 
but they are extremely rare, and contribute marginally to 
the presently known A^BH-A^host relation. In particular, 
the 2 < z < 3 quasars should appear as inactive mas- 
sive galaxies with A4bh~ 10 9 ' 5 Mq in the nearby Uni- 
verse. In order to quantify the occurrence of such objects 
in the Local Universe, we assume the mass function of 
quasar s: $(Xbh ~ 10 9 - 5 M B ,2 < z < 3) « 4 x 10" 9 Gpc" 3 
Mq -1 (|Vestergaard fc Osmerllioogh , and take a volume cor- 
responding to the most distant inactive BH for which a di- 
rect measurement of the mass is availabl^E Under these as- 
sumptions, we expect virtually no objects (0.2 in the whole 
volume) with such high values of V. 

Using the same arguments for targets at intermediate red- 
shift (1.0 < z < 1.5), we expect few (~ 2) objects. Quasars 
at z ~ 1.2 have F values ~ 0.3 dex larger than that at z = 0. 
This offset is close to the one observed for the handful of ob- 
jects populati ng the high-mass end o f the local A^BH-A'fhost 
relation (e.g. iMarconi fc Hundl2003T ). This scenario is thus 
consistent both with the T—z relation of quasars and with 
the observed shape of the local AlBH-A'lhost galaxy relation 
of nearby inactive galaxies. 



5 CONCLUSIONS 

In this paper, we studied the AlBH-A'lhost relations as a 
function of redshift in a sample of 96 quasars from the 
present age to z = 3, i.e., 85 per cent of the Universe age. 
We found that the A^BH-A'lhost ratio increases by a fac- 
tor ~ 7 from z — to z = 3. This trend is not affected 
by significant contributions due to target selection criteria 
and observational biases. Moreover it is independent of the 
quasar luminosity and of the radio loudness. 

1 The BH at the cen tre of the brightest cluste r galaxy in Abell 
1836, at 147 Mpc; see lDalla Bonta et al.l ||2009|'| 



We interpret this trend as an indication that the most 
massive black holes, living their quasar phase at high red- 
shift, keep their high V down to the present age, becoming 
very rare objects in the Local Universe. A fully consistent in- 
terpretation of these results in terms of the common history 
of black holes and galaxies requires further efforts in refin- 
ing the picture sketched here. In particular, two key points 
are yet to be clarified: i- how the mechanisms of quasar 
feedback act onto the host galaxies; ii- which is the role of 
both dry and wet mergers concerning the quasar activity 
and in triggering star formation. Moreover, a better knowl- 
edge of the A^BH-host galaxy relations (improving statistics 
at high masses) will clarify whether very massive, quiescent 
black holes can actually be found in galaxies in the Local 
Universe. 
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